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We derive an exact expression for the time-averaged electromagnetic energy inside a chiral dispersive sphere ir- 
radiated by a plane wave. The dispersion relations correspond to a chiral metamaterial consisting of uncoupled 
single-resonance helical resonators. Using a field decomposition scheme and a general expression for the elec- 
tromagnetic energy density in bi-anisotropic media, we calculate the Lorenz-Mie solution for the internal fields 
in a medium that is simultaneously magnetic and chiral. We also obtain an explicit analytical relation between 
the internal electromagnetic field and the absorption cross-section. This result is applied to demonstrate that 
strong chirality leads to an off-resonance field enhancement within weakly absorbing spheres. 
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1. Introduction 

Metamaterials are artificial structures with engineered 
electromagnetic (EM) response that may exhibit un- 
usual properties such as negative refraction [J, image 
resolution beyond the diffraction limit @, optical mag- 
netism Jl, 0| , EM cloaking [jjj, Q , and slow light propa- 
gation |7| . One possible route to achieve negative refrac- 
tion requires that the real parts of the electric permittiv- 
ity and magnetic permeability are simultaneously nega- 
tive [l|, ■ This condition can be achieved with metama- 
terials composed of periodic structures containing both 
electric and magnetic resonators [ij. Another route for 
negative refraction involves single-resonance chiral meta- 
materials (CMMs), with neither electric permittivity nor 
magnetic permeability negative required Q. Although 
physically different, both routes explore resonant meta- 
materials, which are necessarily absorptive and disper- 
sive. 

The EM energy density in absorptive and dispersive 
media, which is crucial for many applications in photon- 
ics, was firstly derived for nonmagnetic material s [9|| and 
then generalized for composite magnetic media [10( . For 
the particular case of wire-split ring resonator metama- 
terials, there basically exist two methods to calculate the 
EM energy density: the equivalent circuit [ll| and the 
electrodynamic [12[ approaches. These methods were 
shown to be non-equivalent and the apparent inconsis- 
tency between them has been solved in [13]. Specially, 
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this latter result has been recently used to calculate the 
EM energy within a coated magnetic sphere TJ] , an im- 
portant geometry for applications in plasmonics. The 
methodology developed in [13j has been also employed to 
derive the EM energy density stored in single-resonance 
CMMs composed of uncoupled helical resonators flij . 
To the best of our knowledge, the case of the EM en- 
ergy within a dispersive sphere which is simultaneously 
chiral and magnetic, such as the CMM composites, has 
not been treated so far. 



The aim of this paper is to derive an exact ex- 
pression for the time-averaged EM energy inside a 
CMM sphere and its connection to the absorption cross- 
section. Based on the Bohren's decomposition of the 
EM field for optically active spheres [l6], we obtain the 
Lorenz-Mie solution for the internal fields in the mag- 
netic case [13, G3 . We investigate the behavior of inter- 
nal resonances inside the sphere (near-field) and their 
corresponding resonances in the extinction efficiencies 
(far-field). In particular, we derive an explicit expres- 
sion for the absorption efficiency in terms of the internal 
EM fields that reveals that strong chirality leads to an 
off-resonance field enhancement within weakly absorb- 
ing spheres. 

This paper is organized as follows. In Sec. [31 we for- 
mulate the problem of EM scattering by spherical CMMs 
and derive the main analytical results of the paper. In 
Sec. |31 we present the numerical results, whereas Sec. 0] 
is devoted to the conclusions. 
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2. Basic theory 

Let us consider a plane and monochromatic EM wave 
[E(r), H(r)] T exp(— lut), where uj is the angular fre- 
quency and T is the transpose operator. This wave is 
incident to a linear, spatially homogeneous and isotropic 
CMM sphere, with radius a and dispersive optical prop- 
erties [ei(w), Hi(ui), k(ui)], where t\ and \i\ are the elec- 
tric permittivity and magnetic permeability, respec- 
tively, and K is the dimensionless chirality parameter. 
The surrounding medium is assumed to be the vacuum 
(eo,/Xo) and the constitutive relations concerning to the 
sphere material are 



Di = ei(w)Ei +z^M Hl , 

c 

c 



(1) 
(2) 



where c = l/^o/xo) 1 / 2 is the speed of light in vac- 
uum and i = — 1. The refractive indices associ- 
ated with the right-circularly polarized (RCP,+) and 
left-circularly polarized (LCP,— ) waves are rrVy = 
( £ i/ e o) 1 ^ 2 (^i/A t o) 1 ^ 2 ± k. Hence, for a real positive k, it 
is possible to obtain negative refraction for LCP waves 
even if the real parts of ei and /xi are positive Q. 

The constitutive relations in Eqs. (JIJ and ([2]) for 
time-harmonic fields are based on the symmetrized Con- 
don relations as discussed by Silverman [lj| |2(j. Fol- 
lowing [Uj], the constitutive parameters of a single- 
resonance chiral metamaterial are 1211: 



ei(w) = e 



1 - 



UJ 2 — UJ 2 , + iTuj 



1 - 



Flo 2 



U! 2 — UJq + iTlu 

Auj 



k(uj) = - 2 

U1 Z — Uln + 11 UJ 



(3) 

(4) 
(5) 



where uj p and loq are the plasma and resonance frequen- 
cies of the resonators, F is a dimensionless filling factor 
in one unit cell (0 < F < 1) , T is a dissipative coefficient 
and A = ±uj p F 1 / 2 [H, |H|. Both pb\ and k depend on 
the filling factor F; in particular, \k\ 2 oc F. To guar- 
antee the validity of these constitutive parameters, we 
assume that the relevant wavelength within the sphere 
is much longer than the lattice constant associated with 
one unit cell of the uncoupled helical resonators [2l|. 

2. A. Energy density in single-resonance chiral meta- 
material 

For a harmonic wave in a medium [ei(u>), hi(uj), k(uj)], 
described by Eqs. ([l])-©, the time-averaged EM energy 
density is 



<«>, = ~ 



e cff |Ei |^ + McfflH!! 2 + -Re ( Kcff EJ • Hi) 
c 



m 



where the effective energy coefficients are given by [13j, 



e e ff(w) = £o 



/ieff(w) = /io 



UJ* [UJ^ + UJ 



(wg - u 2 y - 

Fuj 2 (3uj 2 



T 2 uj 2 



{uj 2 -uj 2 Y + Y 2 l 



Keff(w) = - 



Auj (Tuj 



2iujI) 



{uj 2 - uj 2 Y + Th 



Under the lossless assumption (r = 
efficients assume the well-known form 



d(ujfii)/duj 

,2 



22] 



(7) 



(8) 



(9) 



these co- 
e c ff = 



> and k 



cir 



d(ujei)/duj > 0, fj, c s 

id(ujK)/duj, with e c ff/i c ff > k 2 s . Therefore, the time- 
averaged EM energy (W)t within the metamaterial 
sphere is obtained from the integration of Eq. ([5]) with 
respect to the volume of the particle: 

(W) t = d<j> d(cos0) / drr 2 (u) t , (10) 
Jo J-l Jo 

where r is the radial variable, and 9 and <p are the polar 
and azimuthal angle variables in a spherical coordinates 
system centered in the sphere, respectively, and (u)t — 
(u) t (r, cos 8, cf)) is given by Eq. (J6)). We emphasize that 
Eqs. ([7])-([n]) are not general, since they are derived for a 
special kind of CMM medium with constitutive relations 
and parameters given in Eqs. ([lj-©. 

2.B. Electromagnetic scattering by an optically ac- 
tive sphere 

The exact solution for the EM wave scattering by an op- 
tically active sphere is developed in [l6[. Here, we gen- 
eralize this solution to dispersive media that are both 
chiral and magnetic, and explicitly calculate the mul- 
tipole moments associated with the internal EM field. 
Following [l(|, the sphere constitutive equations are 

(DBF) 



(DBF) 

Mi 



(Ei 4 
(Hi 



0VxEi) , 
-/3V x Hi) 



(11) 
(12) 



which are known as Drude-Born-Fedorov (DBF) rela- 



tions 



The quantities e^ DBI 



M 



(DBF) 



and 



are 



phcnomcnological coefficients that describe the optical 
activity. Equations (fTTj) and (|12j) are different from the 
corresponding constitutive Eqs. (P) and |2]). However, 
assuming the time-harmonic dependence e~ lu}t , the curl 
macroscopic Maxwell's equations V x Ei = jwBi and 
V x Hi = -jwDi hold for both set of constitutive re- 
lations. Setting a = /3uj and 7 = k/c, one can easily 
demonstrate that the Condon and DBF constitutive re- 
lations are equivalent by the substitutions [H, [24j 

„,2 

__(DBF) 



(DBF) 
PI 



7 

= ei 

Mi 

7 2 

= Mi 

ei 



eifix - Y 



(13) 

(14) 
(15) 
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with inverse relations [24[: ei = < 

(DBF) (DBF) 



^(DBF) i ^ an ^j ^ _ 

i-a 2 4 DBF vr F) - 



Mi 



(DBF) , > 
1 /<■>! Ml 

/£, where £ 



The relations above between the constitutive param- 
eters in DBF approach and Eqs. (P) and © allow us to 
write the Lorenz-Mie quantities in the system of param- 
, /3] and calculate the time-averaged 



eters e 



(DBF) 



(DBF) 

Mi 



EM energy in the system of parameters (ei,/ii,/e). Ac- 
cording to ) we can relate these phenomenological 
coefficients to the complex refractive indices TOl and tor 
through 



P 2 U R k L 



.(DBF) (DBF 

: i Mi 



1/2 






.2 (fc R + k L ) 



(16) 
(17) 



where k q = km q is the wave number in the active 

medium [ep BF ' ) , p,^ BF ^ , /?] , with q is either L or R 
for refractive indices associated with left-circularly po- 
larized (LCP, — ) or right-circularly polarized (RCP,+) 
waves, respectively, and k is the incident wave num- 
ber. Explicitly, one has l/fc R 

;F1 ■ i DBF. .. 

Mi 



l/fc< DBF) + f3 and 



( DBF ) = C jfe (DBF) /t. (DBF) l 1 / 2 



1/feL = l//cP BF) -/3, with k\ 
Since the real parts of g^ DBF - ) and /ip BF ' ) can both as- 
sume negative values, it is convenient to write fc 1 DBF '* = 

w[e (DBF) ]1/2[M (DBF) ]1/2 _ 

The macroscopic Maxwell's equations for harmonic 
waves in DBF approach are (V 2 + /C 2 )(Ei,Hi) T = 
(0,0) T , V x (EkH^ = /C(Ei, Hi) T , V • (Pi,Hi) T = 
(0,0) T , where HU 



nates system (r, 9, 4>) the electric field components: 



E\ r — E n s'm6n(n + 1)tt„ 



^n(pL) r, . , , , n 
2 [Jon Sin (f> + J cn COS 0J 

^(Pr) r • , . 

-an 2 — Son sin (p + g en cos < 

P'r 



(18) 
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PL 
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PR 
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PR 
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-zaei /ij 
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Mi 

(DBF) (DBF) 

-iae 1 p 1 



A^KA = 



To diagonalize /C and simplify the system of equa- 
tions, one can use a linear transformation of the 
EM field 0, [H: (Ei, Hi) T = -4(Q L , Qr,) t , with 

fc n L ? YwhareW 1 « R 
-kn J \ «l 1 

-z[ Ml DBF) /4 DBF) ] 1/2 and a L = -l/a R . The EM left- 
handed and right-handed waves Ql and Qr, respec- 
tively, are the vector basis for the Bohren's decom- 
position scheme and independently satisfy the vector 
Helmholtz equation (V 2 + fc 2 )Q = 0, with V x Q = fcQ 
and V • Q = 0, where k = ki, for Q = Ql and k = — kn 
for Q = Qr. Solving these equations for the trans- 
formed fields Q and returning to the original internal 
EM field (Ei, Hi) T , one obtains in the spherical coordi- 



where E„ = Eoi n (2n + l)/[n(n + 1)], p q — k q r, with 
q is either L or R, and Eq is the amplitude of the in- 
cident EM wave. The corresponding magnetic field Hi 
is obtained from Eqs. (pT8")) - (j2"0)) by replacing (/„, «r<m) 
with (a L f n ,g n ). The radial functions ip n (p) = pj n (p) 
and £ n (p) = p[jn(p) + Wn(p)] are the Riccati-Bessel 
and Riccati-Hankel functions, respectively, where j n and 
y n are the spherical Bessel and Neumann functions. 
The angular functions are 7r„ = (cos 6*)/sin 8 and 
r n = dP^(cos9)/d6, with P^cos^) being the associated 
Legendre function of first order. Expansions of the inci- 
dent and scattered EM fields in terms of vector spherical 
harmonics can be found in 1251. 



Using Bohren's notation [lfj, the Lorenz-Mie coeffi- 
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cients for the splitted internal fields are: 



Jon 
fen 
9on 



irrihWi 



(R) 



(R) 



V n ] = 1pn{mqX)£' n (x) - fh^ n {x)lp' n (m q x) , 

W n q) = mi/j n (m q x)£' n (x) - £ n (x)ip' n (m q x) , 



(21) 
(22) 
(23) 

(24) 

(25) 
(26) 



where x — ka is the size parameter and the effective 
impedance index fh is 



m 



(DBF) 



"I 1/2 



(DBF) 



Mo 



(DBF) 

Mi 



f 



I 

m L 



with m = m/i^ DBF ' ) //zo being the effective re- 
fractive index associated with the chiral medium 
[e[ DBF \/z^ DBF \/3]. In particular, if the medium is non- 
magnetic [^^ DBF ' 1 = /in], one has m = fh (25[. The inter- 
nal coefficients g on and g cn have dimensions of «l and, 
thereby, it is expected to appear in the dimensionless 
calculated quantities the products a^g on and a^g cn . 



2.C. Electromagnetic energy within chiral spheres 

To calculate the time-averaged EM energy within a 
spherical particle, we follow the same procedure of [bj . 
I2a - l28j |. We define the partial contributions to the aver- 
age EM energy (W)t, since this quantity can be written 
as a sum of the electric, magnetic and magnetoelectric 
coupling terms: 



(W) t = (W E )t + (W„) t + (W EH )t 



(27) 



where (W E )t = e eS J v d 3 r|Ex | 2 /4, (W H )t = 
Meff/yd^lH!! 2 ^ and (W EH ) t = J y d 3 rRe( CTcff Hx • 
E*)/2c, with the region of integration V being the 
volume of the chiral sphere. For a sphere with ra- 
dius a having the same optical properties as the 
surrounding medium (en, up), the average EM energy is 
(Wo) t = 27rah \E \ 2 /3 

In particular, consider the radial and angular contri- 
butions to the electric energy: (W E ) t = (W Er ) t + (W E g + 



W E ^,)t- Explicitly, the radial contribution is 
(W Er ) t = ^^ dtj d(cos6)J dr r 2 \E lr 



T^(Wo)t— $>(n + l)(2n + l) 



4a 3 



(l/on| 2 + \U 



M 2 (l 



gon 



\9en\ 



\Upl)\ 2 

|fc L | 2 

\jn( P n)\ 2 



2Re 



a R {fang on + fcng 





in(PL)in(P R ) 



(28) 



where we have used (2n + I) f_ t d(cos 0)ir n iT n i sin 2 9 = 
2n(n + l)5 n ,n> to simplify this expression [IHH^]- Anal- 
ogously, the angular contribution to the electric energy 
is (Web + W E4> ) t = e off / v d 3 r(\Eie\ 2 + \E^\ 2 )/4, which 
leads 



(Weo + W E *)t 



(W )t— V(2n+1) 

t=r, ^ ' 



4a 3 



n=l 




|2,|/ |2\ \4>n(PL)\ + K(Pl)T 
\\Jon\ + |/e: 

l«R| 2 (\9on\ 2 + |3e„| 2 ) 



|fc L | 2 



\kn\ 2 



2Re 



a R (/onSon + fen9en) 



^n(PL)^ ra (p R ) ~ 1pn(pL)4>n(PB) 



(29) 



where we have used the following relations |27j |: 
(2n + f) f\ d(cos6)(TT n n n , + T n T n /) = 2n 2 (n + l) 2 5 n , n > 

and j\ d(cos 6)(n n T n > + T n n n ,) = 0. 

From recurrence relations involving spherical Bessel 
functions, (2n+l)j n (p) = p[j n -i(p)+jn+i(p)} and (2n+ 
!)in(p) = n 3n-i(p) - (n + l)j n+1 (p) [27], we obtain 



(2n+l) 
PAP* B 



[n(n + l)jn{PA)3n{p*B) + C(Pa)C(/°b)] 



nj n+ l(p A )in+l(PB) + (n + l)jn-l(PA)in-l(PB) 



(30) 



where A and B can be either L or R. Note that the 
left-hand side of Eq. (|30|) appears when we calculate the 
average electric energy (W E ) t as the sum of Eqs. (|2"g)l 
and (|2"9"|) . Indeed, the right-hand side of Eq. (I5TJ1) is used 
to simplify the calculations, once we have the analytical 
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solution for the integral [271 . |2S 
dr r 2 j n (p A )j n (p* B ) 



3 [VBjn{VA)3n{yB) ~ Z/Ajn (?/b) (2/a)] 



Vl Vb 2 



(31) 



with ua — rriAX, ub = tubx, and we assume m\ ^±mg. 
For m A = ±m£, the L'Hospital's rule provides [27| | 



lim 

tua— >±mp 



dr r 2 j n (m A kr)j n {m* B kr) 



,3,(-l±l) 



[j'n(yA) -Jn-l(2/A)jn+l(2/A)] 



(32) 

Therefore, defining the dimensionless functions 

2~i AB) = 4 / dr r 2 j„(m A fcr)j n (m* Ax) , (33) 



4 A ± B) = nZ#? + (n + l)li AB) ± (2n + 1)1^ (34) 



r(AB) 



(AB) 



^ (±) =E (l/on| 2 + l/e»| 2 )4 L + L) 
n=l 

+ |aR.| 2 (l5on| 2 + Iffcnl 2 ) •^l R +' ) 



"R (/onffon + feng* en ) FT- } , (35) 



±2Rc 



and summing (W E )t = (W E r)t + (W Eg + W E< j>) t , from 
Eqs. (gSJ) and and using Eq. ([201), we finally obtain 



(Wg)t = 3£eff 
(W )t 4 e B ' 



(36) 



where S E — S^- ' is a sum associated with the volume 
integral of |Ei| 2 . Analogously, from the calculations 
of (W ff )t = (W Hr ) t + (W„ e + W H $)t and (W EH )t = 
(W EHr )t + (W EHe + W B H<t>)t, we obtain the magnetic 
and magnctoclcctric coupling terms: 



{W H )t 3 2 Mefif 

/T „ . = t m 6h, 



(Mj = 3 
<VF ) t 2 

where 5// = and 



Im (mK e eS EH ) 



(37) 
(38) 



5eH — E <j (\fon\ 2 + |/on| 2 j 
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2 / ,2,| ,2\ r (RR) 
-|«r| (Jffon| +|.9cn| J-^„,+ 

«R (/onffon + fcngDJ 7 ^ > (39) 



-2dm 



are related to the volume integrals of Hi | 2 and (E| -Hi), 
respectively. In Eq. (|39l) . we have considered the fact 



that (^o/eo) 1 ^ 2 = ima^. Equations (l36l) ~ (j38)) are the 
main analytical results of this study. Together with the 
constitutive relations and parameters associated with a 
single-resonance CMM, Eqs. ©-© and Eqs. 
and the transformation relations from this parameters 
system to the DBF approach, Eqs. (p~5|) (fX5T) . we can 
readily calculate the average EM energy (W)t within a 
CMM sphere. In particular, one can obtain the time- 
averaged power loss {P)t from the internal energy {W)t 
by means of the following substitutions in Eqs. ((361) 
(EED 0: 



e c ff 



u)e"((J) = e 



Tuj^uj 



2, ,2 



— > IUJK (OJ) = 

2 v ; 



TFuj 4 



-iTAuj 3 



(w 2 - oj 2 Y + r 2 w 2 



(40) 
(41) 
(42) 



where (e", ///, k") — Im(ei, fii, k), i.e., the imaginary 
parts of the constitutive quantities. Also, the scatter- 
ing, extinction and absorption efficiencies are calculated, 
respectively, as follows [25| : 



Qsca,± 

<9ext,± 
<5abs,± 



n— 1 

±2Im [(o n 

_ oo 
n— 1 

Qcxt,± — Q sca,± 



l){l«n| 2 - 

K)c* n ]} 
1) Re (a„ 



|6„| 2 + 2|c„| 2 

(43) 

- b n ± 2ic n ) , (44) 
(45) 



where the signal "+" is chosen for RCP waves and "— " 
for LCP ones, and the scattering Lorenz-Mie type coef- 
ficients a n , b n , c n and d n are: 



WW 



d„ 

Ail) 



mip n (m q x)ip' n ( x ) ~ ^n{xW n {m q x) 
i/) n (m q x)i/)' n (x) - mip n (x)ip' n (m q x) 



(46) 
(47) 

(48) 

(49) 
(50) 
(51) 



The mean absorption efficiency Q a b s due to both RCP 
and LCP waves is given by Q abs = (<3abs,L + Qabs,R.)/2. 
Following [13, the power loss {P) t is proportional 
to Qabs: and, thereby, depends on S E , Sh and S E h 
calculated above. We have analytically obtained that 
the exact relation between the internal fields and the 
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absorption efficiency is 



^abs ^ 



= 2x 



(W E ) t 



e e s (W )t 



< {W H ) t 3 p 

Meff (Wo)i 2 



£0 



Iml 2 — S H + 2Re (mK"S EH ) 
Mo 



,(52) 



where x = fca is the size parameter. For weakly absorb- 
ing spheres with positive refractive index m = ml + im" 
(m" < to'), it is well-known that (W) t /(W )t ~ 
{8/3)m'x{Q ahs /m") 0, [27j]. Therefore, in the weak 
absorption regime, one can approximately calculate the 
behavior of the internal energy (W)t from Q a bs- 

3. Numerical calculations 

Here we present numerical calculations obtained from a 
computer code written for the free software Scilab 5.3.3. 
According to Zhao et al. [2l[ , we choose realistic param- 
eters for the CMM sphere in the frequency range of Ter- 
ahertz: ujq = ujp = 2 THz and T = 0.05loq. The validity 
of the effective constitutive quantities (ei, k) is guar- 
anteed for size parameters x = ka < 1. For io ~ 10 12 Hz, 
ka < 1 is satisfied for a ~ 10~ 4 m. 



at 



2.0 


(a) 








F=0.\ 






F= 0.5 


1.6 




F=0.9 - 


1.2 






0.8 






0.4 







0.2 



0.4 0.6 

/(THz) 



0.8 



-(b) 






F=0.1 




1 F=0.S 




1 F=0.9 







0.2 



0.4 0.6 

/(THz) 



Fig. 1. EM scattering by a CMM sphere with radius a — 
10 -4 m and fo = ljo/(2tv) ~ 0.32 THz (w p = wo) as a func- 
tion of the frequency, (a) The real part of the LCP refractive 
index ttil for three chirality parameters k(lj) with filling fac- 
tors F = 0.1, 0.5, 0.9. (b) The corresponding LCP extinction 
efficiency Q e xt,L- There is no negative refractive indices for 
these parameters. 
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Fig. 2. EM scattering by a CMM sphere with radius a = 
10 -4 m and fo = u)q/(2-k) ~ 0.32 THz as a function of the 
frequency, (a) The real part of the RCP refractive index 
tub. for three chirality parameters k(u>) with filling factors 
F = 0.1, 0.5, 0.9. (b) The corresponding RCP extinction 
efficiency Qext.R- The negative refractive indices for RCP 
waves occur for /> Jo- 



in Figs. [TJ and O we show the extinction efficiencies 
and the associated refractive indices for LCP and RCP 
waves, respectively, corresponding to EM scattering by a 
CMM sphere with effective parameters given by Eqs ([3])- 
([5]). Three values of the filling parameter F are consid- 
ered here: F = 0.1, 0.5, 0.9. Note that the real part of 
TOl in Fig. HJa) does not assume negative values in this 
frequency range and decreases with the filling parameter 
F for / < fo ~ 0.32 THz, and increases for frequencies 
right above f (f = 0.4 THz to 0.5 THz). As it can 
be noted in Fig. Hlb), Qcxt,L decreases with increasing 
the chirality parameter \k\ cx F 1 / 2 . The opposite is ob- 
served in Fig. [5] For frequencies above fo « 3.2 THz, 
the real part of tor is negative and |torJ increases with 
the chirality parameter k, Fig. Ufa). This leads to an 
enhancement in Q C xt,R, Fig. d{b), where the resonance 
peak for Re(TOR.) < is displaced to high frequencies. 

An interesting result is obtained when we calculate 
the internal energy, Fig. [3ja). Note that, although the 
mean extinction efficiency Q ext = (Qext.L + <3ext,H.)/2 
in Fig. Ufb) is increasing with the chirality parameter 
k for / > fo and Re (tor) < 0, the internal energy is 
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Fig. 3. EM scattering by a CMM sphere with radius a — 
10 -4 m and fo — u)q/{2-k) ~ 0.32 THz as a function of the 
frequency, (a) Average energy inside the CMM sphere for 
filling parameters F — 0.1, 0.5, 0.9. (b) The corresponding 
mean extinction efficiency Q cx t- The decrease in (W)t for 
/ > fo coincides with the increase of Qcxt- 



decreasing simultaneously. For a sufficiently strong chi- 
rality parameter, the LCP waves tend to be suppressed, 
as can be seen in Fig. [Ha). Correspondingly, strong chi- 
rality favors the reflection of RCP waves on the interface 
between the CMM sphere and the surrounding medium. 
Since Q oxt increases as k increases for Re(mR.) < 0, we 
conclude that incident waves are strongly reflected by 
the surface of the sphere and, for this reason, do not 
contribute to the internal energy. 

In Fig. 21 we use the same parameters as above, as- 
suming now F — 0.8 ("strong" chirality) as a function 
of the absorption coefficient V = Gujq, where G = 0.05, 
0.005, 0.0005. Here, we notice some regions where 
weak absorption (Q a bs ~ 0) leads to weak extinction 
(Qext ~ 0), and increases the amount of energy stored 
in the CMM sphere. This effect occurs above fo, in 
the range 0.37 THz< / < 0.40 THz, see Figs, gja) and 
Bfb). For / fts 0.37 THz, Q a b s and Q ext almost vanish 
for G = 0.0005, but the value of (W) t in Fig. gfc) is 
still large. In addition, there are resonance peaks in the 
internal energy around fo that do not correspond to a 
peak in Qcxt, such as the one at / w 0.28 THz. This 
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Fig. 4. CMM sphere with radius a = 10 -4 m, filling factor 
F — 0.8 and ljq — u p — 2 THz as a function of frequency and 
the absorption parameter G = 0.05, 0.005, 0.0005. (a) The 
mean absorption efficiency Q a b s - (b) The mean extinction 
efficiency Q cx t. (c) The internal energy (W)t. For / > /o ~ 
0.32 THz, a decrease in Q a bs and Qext do not necessarily 
coincide with a decrease in {W)t- 



behavior is indicated by a solid arrow in Figs. Bib) and 
@Jc) and suggests the existence of off-resonance field en- 
hancement within the sphere. These resonances are re- 
lated to interferences between two different Lorenz-Mie 
type coefficients, such as f n g* in (W) t , Eq. (|3"5")) . and 
a„c* in Q sca , Eq^ (|43l) . and are explained in terms of 
the Fano effect I31II. This effect also occurs for coated 
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spheres containing dispersive Drude-like materials with 
low absorption [321 ] - 

4. Conclusions 

We have calculated the time-averaged energy inside a 
CMM sphere using a set of dispersive realistic parame- 
ters provided by an effective medium theory [2l| . Us- 
ing the decomposition of the EM field in LCP and 
RCP waves [16[, we have obtained the internal fields 
for the magnetic case and the respective Lorenz-Mie co- 
efficients in DBF approach. Since the energy density in 
an isotropic chiral media is commonly calculated from 
the Condon constitutive relations [l9| , a transformation 
of parameters has been applied to obtain the exact ex- 
pression for the internal energy. We have expressed the 
absorption efficiency in terms of the internal energy and 
the constitutive parameters. Although the internal en- 
ergy decreases with the chirality parameter in the region 
of negative refraction for RCP waves, the mean extinc- 
tion efficiency increases as the chirality parameter in- 
creases. Our results reveal that, in the small-particle 
limit, off-resonance field enhancement occurs for weakly 
absorbing spheres with strong chirality. We show that 
this effect is due to field interferences between different 
multipolar moments, which lead to Fano resonances in 
electromagnetic scattering. We hope that our findings 
could be exploited in applications involving dispersive 
CMM scatterers and Fanoshells. 
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